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1. Introduction 

The mechanism of dynamical chiral symmetry breaking and confinement is one of the most 
fundamental problem of hadron physics. The propagator of dynamical quarks in the infrared region 
provides information on dynamical chiral symmetry breaking and confinement. In the previous 
paper [j|], we measured gluon propagators and ghost propagators of unquenched gauge configu- 
rations obtained with quark actions of Wilson fermions (JLQCD/CP-PACS) and those of Kogut- 
Susskind(KS) fermions (MILC) [§, ||] in Landau gauge and observed that the configurations of the 
KS fermion are closer to the chiral limit than those of Wilson fermions. 

In the analysis of running coupling obtained from the gluon propagator and the ghost propa- 
gator, with use of the operator product expansion of the Green function, we observed possible con- 
tribution of the quark condensates and A 2 condensates in the configurations of the KS fermion[jl]]. 
The quark propagator of quenched KS fermion was already measured in [Q], and possible contri- 
bution of these condensates are reported. Unquenched KS fermion propagator of 20 3 x 64 lattice 
(MILC C ) was measured in [^J, but to distinguish the gluon condensates and the quark condensates, 
it is desirable to measure the quark propagator of larger lattice (MILC f) and to compare with data 
of MILC f . We measure quark propagator of gauge configuration of 1) MILC C 20 x 64, j8 = 6.76 
and 6.83 and 2) MILC/ 28 3 x 96, j3 = 7.09 and 7.11, using the Staple+Naik action|]]. 

2. The running coupling in MOM scheme 

The running coupling in MOM scheme is given with use of the vertex renormalization factor 
Z\ as 

a s {q) = a R (n 2 )Z R (^,n 2 )G R (q 2 ,n 2 ) 2 = ^^\z{ q 2 ,p)G(q 2 (2.1) 

where Z and G are the gluon and the ghost dressing function, respectively, and Z\ is the vertex 
renormalization factor. 
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Figure 1: The logG(qa) as a function of 
log qa (a is the lattice spacing of MILC/) of 
MILC/ P imp = 7.09(stars) and MILQ j5 imp = 
6.76(triangles). 
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Figure 2: The running coupling a s (q) as a 
function of log 10 g(GeV) of MILC/ /3, mp = 
7.09(stars) and 7. 11 (diamonds). 
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Our data suggests that the gluon propagator is infrared finite as in Dyson-Schwinger equation[Q] 
The running coupling in the infrared is suppressed as shown in Figure 2,but the main origin is the 
suppression of the ghost propagator in the infrared. 



We parametrize the difference of the lattice data and the pQCD 4-loop result[]10|, |11|] in the 
lGeV< q < 6GeV region in the form, with a minor correction term d as 

Aa s (q) = a sM ,{q) - a SiPert {q) = — + — +d, (2.2) 

where the A 2 condensates gives c\ and the gluon condensates and/or quark condensates gives c 2 . 
Although statistics is not large, running coupling of CP-PACS suggests c\ ~ 2GeV. The MILC data 
suggests c\ ~ 4GeV and c 2 ~ — 2GeV. There is an analytical calculation that suggests correlation 



between the A 2 condensates and the the horizon function parameter! 12]. 

3. The quark wave function renormalization 

We renormalize the quark field as Ybare = VZiWr> an d define the colorless vector current 
vertex [[7]] 

T^q,p) =S- l (q)G^q,p)S- i (q + p) (3.1) 

where 

G^(p,q) = J ^V >3,+i '^(?(y)?Wy M 9W?(0)} (3.2) 

and S(q) is the quark propagator. 

The vertex of the vector current with p = is written as 

^ti(q) = 8 a ,b{gi(q 2 )Yn+ig2(q 2 )Pn + g^q^q^q + ig4(q 2 )[Y^,4}} (3-3) 
The Ward identity implies 

Z f^T ll {q) = -i^S-\q) (3.4) 

where Zy 0M g\ (q 2 ) = Z v (q 2 ) and Zy OM = 1 when there is no lattice artefact. 
The running coupling g of the ghost, anti-ghost, gluon coupling 

g(q) =Z^zl /2 (ll 2 ,q 2 )Z 3 (ll 2 ,q 2 ) (3.5) 

and that of quark, gluon coupling 

g(q) =Zi l zl /2 (ll 2 ,q 2 )Z 2 (ll 2 ,q 2 ) (3.6) 



Z 2 (H 2 ,H 2 ) = 1 and Z 3 (n 2 ,n 2 ) = 1, and so Z^ l Z^ /2 Z 3 = Z^ l zl /2 Z 2 implies Z\ = Z y . On the 



are identical due to the Slavnov-Taylor identity. At the renormalization point q = fi, we fix 

and Z 3 (n 2 ,n 2 ) = 
other hand Zy(q 2 ) = Zy 0M g\ (q 2 ) where for 16 tastes 

^) = 4^ Tr ^P = °Xtt Jr)] (3-7) 



When Zy 0M = 1, gi(n 2 ) is identical to Z\ which is defined by the renormalization of the running 
coupling on the lattice defined at \i ~ 6GeV and summarized in Table [[]. 
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Pi/K sea i 


fh/Ksea2 


average 


configurations 


CP-PACS 


1.07(8) 


1.21(10) 


1.14 


W 2 = 0.1357,0.1382 


MILC C 


1.49(11) 


1.43(10) 


1.46 


Pi )2 = 6.83,6.76 


MILC/ 


1.37(9) 


1.41(12) 


1.40 


01,2=7.11,7.09 



Table 1: The 1 jl\ factor of the unquenched SU(3). 
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Figure 3: The renormalization factor Z%(q) of 
MILC/- with bear mass mo = 13.6MeV(stars), 
27.2MeV(diamonds) and 68.0MeV(triangles). 
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Figure 4: Same as Fig. 3 but quark of 
MILQ with bear mass m = 1 1 .5Me V(stars), 
65.7MeV(diamonds) and 82.2MeV(triangles). 
The last two almost overlap. 



The MILC configurations are produced by 1 _ 7 5 
the Asqtad action with the bare masses 13.6MeV, i . 5 

27.2MeV and 68.0MeV in the MILC /, and 1 1 .5MeV, 1.25 
65.7MeV and 82.2MeV in the MILC C . The gluon 
wave function renormalization Z3 of the Asqtad 
action is renormalized by 1/mq, where uq is the 
fourth root of the plaquette value. The tadpole 
renormalization on quark wave function is uq and 
we normalize S(q) by multiplying the average 
of w Zi i.e. 1/1.36 for MILC/ and 1/1.38 for 
MILC C . The quark propagators 7si{q) after this 
renormalization are infrared suppressed as shown 
in Figures 3 and 4. The apparent difference in the 
formulae of [^|] and our work are only in the ex- 
pression and in fact they are equivalent and the 
results agree with [g|. 
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Figure 5: The quark wave function renormaliza- 
tion factor Z v (q) of MILC/- with bear mass mo = 
27.2MeV. Dashed line is the pQCD result and solid 
line is the result including the A 2 ^ condensates. The 
points below 0.6GeV are g\(q 2 ) approximated by 
the running coupling obtained by the ghost anti- 
ghost gluon vertex. 



Using the pQCD result of the inverse quark 
propagator[|l(], 11] and (A 2 ) and C2 (the contribution from the mixed condensate (qAq) [Jsj] ) as fitting 



parameters, we calculate 



Z ¥ (q) 



Ml) 



Z\ff 



(l 2 ) + 



a(ji) 

a(q) 



-ro+r A 2 



(A 2 (M)) 



pert 



4(JV?-1) 



(M 2 ) + 



C2 



(3.8) 
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where a (q) are data calculated in the MOM scheme using the same MILC /■ gauge configurations [jlj]. 
We estimate that ci is small. In Figure 5, we show a fit of the MILC/ mo = 13.6MeV data, by tak- 
ing the renormalization point at \i ~ 3.8GeV with use of (A 2 (jit)) ~ 1.6(3)GeV 2 and ca = 0. These 
parameters are consistent with [Q]. 

In q > \.5GeV region, dynamical mass of a quark in pQCD is expressed as[ffl] 



M(q) 



47C 2 d M (qq),[log(q 2 /A 2 QCD )] d ^ m( M 2 ) [log^/A 2 ^)] 



3q 2 [log(^ 2 /A 2 QCD )} 



+ ■ 



[log(^/A 2 CD )]^ 



(3-9) 



where dj^ 



12 



33-2/V, 



. The second term is the contribution of the massive quark. In the analysis of 



the lattice data, we observe that the quark condensates — (qq)^ and Aqcd roughly satisfy — (qq)^ 
(0.70A ecZ )) 3 [|], with A QCD = 0.69GeV. 

For the global fit of M(q), we try the phenomenological monopole typejl3|] 

cA 3 

M{q) = - T - T ^ + m Q (3.10) 



q 2 +A 2 



where mo is the bare quark mass. 

Figures 6 and 7 show the mass function M(0) of MILC/ and MILQ , respectively. 




Figure 6: The dynamical mass of the MILC/ 
quark with bear mass mo = 13.6MeV(stars), 
27.2MeV(diamonds) and 68.0MeV(triangles) 
and the phenomenological fits. 



Figure 7: The dynamical mass of the MILQ 
quark with bear mass mo = 1 1 .5MeV(stars), 
65.7MeV(diamonds) and 82.2MeV(triangles) 
and the phenomenological fits. 



We observe that the product cA becomes larger as the bare quark mass becomes heavy and 
it depends on /3 in the case of MILC/ but not in the case of MILC C . In the case of MILC/ mo = 
13.6MeV, the lowest three momentum points of M{q) are systematically smaller than the other 
points. Ignoring these points we find cA of j8 = 7.09 is smaller than that of 7. 1 1 and that of MILC C 
as shown in Figure |8[ In the chiral limit mo — > 0, we obtain M(0) = 0.35 ~ 0.37GeV, which is 
larger than ^ and that of the Wilson fermion[[l3|] by about 20%. 



4. Discussion and conclusion 



We measured running coupling of unquenched Wilson fermion and KS fermion and the quark 
wave function renormalization factor and mass function of the KS fermion. The renormalization 
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factor Zy obtained as I/Z2 is infrared finite. The Kugo-Ojima confinement criterion favours infrared 
vanishing of Zy and gi(q 2 ) approximated by the running coupling suggests this behavior, but it 
could be a lattice artefact. With infrared finite Z2, infrared vanishing of Z\y is necessary for the 
confinement criterion to be satisfied. 
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